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ABSTRACT 

We present high resolution simulations of a multiple merger of three disk galax- 
ies including the evolution of magnetic fields performed with the JV-body/SPH code 
Gadget. For the first time, we embed the galaxies in a magnetized, low-density 
medium, thus modeling an ambient IGM. The simulations include radiative cooling 
and a model for star formation and supernova feedback. Magnetohydrodynamics is 
followed using the SPH method. The progenitor disks have initial magnetic seed fields 
in the range of 1(T 9 to 1(T 6 G and the IGM has initial fields of 10" 12 to 1(T 9 G. 
The simulations are compared to a run excluding magnetic fields. We show that the 
propagation of interaction-driven shocks depends significantly on the initial magnetic 
field strength. The shocks propagate faster in simulations with stronger initial field, 
suggesting that the shocks are supported by magnetic pressure. The Mach numbers 
of the shocks range from approximately M = 1.5 for the non-magnetized case up to 
M = 6 for the highest initial magnetization, resulting in higher temperatures of the 
shock heated IGM gas. The magnetic field in the system saturates rapidly after the 
mergers at ~ 1CU 6 G within the galaxies and ~ 1CU 8 G in the IGM independent of 
the initial value. These field strengths agree with observed values and correspond to 
the equipartition value of the magnetic pressure with the turbulent pressure in the 
system. We also present synthetic radio and polarization maps for different phases of 
the evolution showing that shocks driven by the interaction produce a high amount of 
polarized emission. These idealized simulations indicate that magnetic fields play an 
important role for the hydrodynamics of the IGM during galactic interactions. We also 
show that even weak seed fields arc efficiently strengthened during multiple galactic 
mergers. This interaction driven amplification might have been a key process for the 
magnetization of the Universe. 
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— galaxies: magnetic fields — galaxies: kinematics and dynamics 



1 INTRODUCTION 



Radio observations have revealed that most late type galax- 
ies in the local Universe - isolated grand design spirals, irreg- 
ulars and dwarf g a laxies - are permeated by magnetic fields 
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strengths in all of these objects do not vary by more than 
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one order of mag nitude from a few /iG in dwarfs (e.g. 
IChvzv et al] 120031 ) to 30 /^ G in the star-forming regions of 
grand-design spiral galaxies |Fletcher et alj|2004l ). Magnetic 
fields have also been observed at redshifts up to z m 2 in 
damped Ly-a systems. These systems, which might be in- 
terpreted as large prog enitors of present-day galaxies (e.g. 
I Wolfe et al.| [l995. 2005), seem to host magnetic fields of sim- 
ilar strength as local late-type galaxies (e.g. iBernet et all 
120081 and references therein) . 

These observations invite the question about the origin 
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and the evolution of the magnetic fields in th e early universe. 
Differ ent scenarios have been suggested: iLesch fc Chibal 
l| 19951 ) have shown analytically that strong magnetic fields in 
high redshift objects can be explained by the combined ac- 
tion of an evolving protogalactic fluctuation and electrody- 
namic pro cesses providing magn etic seed fields (i.e. battery 
processes) . I Wang fc Abel (|2009h performed numerical simu- 
lations of the formation of disc galaxies within an collapsing 
halo imposing a uniform initial mag netic field of 1(T 9 G. The 
initial field grew by three orders of magnitude within approx- 
imately 500 Myr of the evolution. The amplification might 
be due to the combined effects of magnetic field compression 
during the collapse and amplification of the uniform initial 
field b y differential rotation as studied also in lKotarba et all 
(|2009h . These studies indicate, that the growth of magnetic 
fields might be a natural part of the galaxy formation pro- 
cess. 

A key ingredient in galaxy formation studies, however, 
is the consideration of galaxy interactions. Within the stan- 
dard cold dark matter (CDM) models present-day galaxies 
have undergone several major and minor mergers at ear- 
lier epochs of the universe, and thereaf ter continued accret- 
ing gas and smaller g alactic subunits (I White fc Rees|[l978l ; 
IWhite fc F rcnk 199lj) . Interactions of galaxies change their 
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tential is varying rapidly during the interaction. Since 
the gas component is dissipative and most sensitive to 
variations of the gravitational potential, it is strongly 
affected by the interaction and driven towards the 
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estimate the amount of free energy during an interaction 
of two galactic subunits to be proportional to their relative 
velocity squared, i.e. Ef TCC ~ v^ el . Obviously some of this en- 
ergy released during the interaction is converted into ther- 
mal energy of hot gas. High energy particles also carry away 
some of the energy. However, it is reasonable to assume that 
at least some of this energy is converted into magnetic field 
energy during the compression of gas and the formation of 
tidal structures. As the amount of Ef Tee can be very large 
during a major merger, the amount of energy converted into 
magnetic energy might be significant. Moreover, gas which 
is heated by the interaction and driven into the IGM should 
carry magnetic energy out of the galactic units, thus mag- 
netizing the IGM. This process shoul d be similar to what 
was found bv lDubois fc Tevssier! (|2009l ) in their study of the 
formation of dwarf galaxies including magnetic fields and 
galactic winds. 

So far, simulations of interactions and mergers of 
disk galaxies have been mainly investigated with re- 
spect to changes in stellar dynamics, gas flows, star for- 
mation (SF ) and the formation of central supermassi ve 
black holes jDi Matteo et alj|2005l; ISpringel et al.l l2005al lbl; 



Since gas and magnetic field are tightly coupled, the mag- 
netic field traces the gas motion and will be strengthened 
by shocks and gas inflow. A perfect example for the strong 
coupling of gas and magnetic fields is the interacting sys- 
tem NGC4038/4039 (th e "Antennae gal axi es") which has 
been r ecently simulated (|Karl et al.ll2010h bv lKotarba et all 
|20ld ) including magnetic fields. It was shown that even an 
initial magnetic field as small as 10 -9 G grows significantly 
during the interaction of two equal-mass spiral galaxies. The 
magnetic field strength thereby saturates at a value of ~ 10 
uG, in good agreement with observations (|Chvzv fc Beckl 
l2004h . This saturation value was reached independently of 
the initial magnetic field stren gth in the range of Bp = 10~ 9 
- 10" 4 G. We emphasize that iKotarba etall (|2010h found 
the saturation level to correspond to near equipartition be- 
tween magnetic and turbulent gas pressure, which is in good 
agreement wi th theoretical consider ations of the turbulent 
dynamo (e.g. lArshakian et al. I l2009l and references therein, 
see also section f4.3.3p . Furthermore. IKotarba et al.l (|2010l ) 
provided synthetic radio maps calculated at time of best 
match between the simulated gas and stellar distributions 
and observations. These synthetic radio maps are in con- 
vincing morphological agreement with synchrotron observa- 
tions of the Antennae system and the underlying numerics 
of the applied JV-body/SPH code Gadget showed to be ca- 
pable of following the evolution of magnetic fields in a highly 
nonlinear environment. 

High resolution simulations of the formation of in- 
dividual galaxies in a full cosmolog i cal context (see e.g. 
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the dramatic impact of mergers on the gas flows will directly 
affect the dynamics of the magnetic field of the systems. 



ing the processes leading to the magnetization of the Uni- 
verse. Although it would be worthwile to consider cosmo- 
logical studies in the long run, the simulations of three in- 
teracting galaxies presented in this paper are an additional 
further step towards a more complete scenari o and a natural 
extens ion of the previous study presented in lKotarba et all 
(|2010l ). With this study we show that the magnetic field 
growth accompanying a galactic interaction and its satura- 
tion at the equipartition level between magnetic and turbu- 
lent pressure holds also for a more general setup of interact- 
ing galaxies including IGM gas. 

The paper is organized as follows: We briefly describe 
our numerical method in section [2] In section [3] we present 
a detailed description of the setup of the three colliding 
galaxies. The temporal evolution of the simulated systems 
and particulary the magnetic fields is described in section 
U In section [S] we present synthetic radio and RM maps 
of our simulated system. Finally, we summarize our results 
and conclude in section [6] 



2 NUMERICAL METHOD 

The simulations presented h ere were perfo rmed with the 
N-body/SPH-code Gadget (|Springelll2005l ). Gravitational 
interactions between th e particles are evalua ted with a hi- 
erarchical tree method (|Barnes fc Hutj["l986T ). The dynam- 
ics of the Lagrangian fluid elements are followed using 
a SPH formulation which conserv es both energy and en- 
tropy (|Springel fc Hernquistl 12002) including the evolution 
of magnetic fields which was implemented and tested by 
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iDolag fc Stasvszvnl (|2009r i. The code has already been used 
to investigate t he evolution of mag netic fields in isolated 
spiral galaxies (jKotarba et alJ l2009h and during the colli- 
sion of two equal m ass spiral galaxies re sembling the nearby 
Antennae galaxies (|Kotarba et aljfefjich . 

The simulations presented in this paper have been per- 
formed using the standard (direct) magnetic field implemen- 
tation. We apply the Lorentz force and artificial magnetic 
dissipation with an artificial magnetic dissipation constant 
of as = 0.5. Artificial magnetic dissipation is included in 
order to reduce numerical errors arising from the SPH ap- 
proximati on rather than to capture any physical dissipation 
correctly (|Dolag fc Stasvszvnl |2009| ). Thereby, only a small 
fraction of the magnetic field within the volume defined by 
the SPH particle is allowed to leave this volume within the 
local dynamical time (defined via the local signal velocity). 
This limitation ensures that even strong shocks are well cap- 
tured. Yet, of course, the artificial dissipation also leads to 
an effective slow diffusion of the magnetic field. This can be 
seen in the simulatio ns of isolated disk galaxies presented in 
iKotarba et"afl (|2010T ) (Fig. 3 in their paper, the yellow and 
red lines). However, the value of as is chosen such that on 
the one hand the numerical errors are efficiently reduced, 
and on the other hand, the magnetic diffusion is prefer- 
ably low. Hence, the numerical dissipation does not reflect 
the true physical dissipation. Physical magnetic dissipation 
arises either due to electric conductivity, in which case it 
is very small and thus negligible, or due to turbulent diffu- 
sion, which can be significantly higher. However, turbulent 
diffusion would have to be modeled within the simulations, 
because it reflects processes on sub-resolution scales. We do 
not include any sub-resolution turbulent diffusion model in 
our simulations. 

In addition to lKotarba et all (|20I Oh . we now also apply 
the subtraction of the effect of numerical magnetic diver- 
gence (the "div ergence force") in the momentum equation 
as su ggested by lB0rve et al.l (|200l[ ) (see IDolag fc Stasvszvnl 
120091 for more detail s). However, as a r efinem ent of the 
method presented in IDolag fc Stasvszvnl <|2009h . we define 
a threshold for the divergence force subtraction: Whenever 
the correction becomes larger than the half of the current 
Lorentz force, it is limited to that level. In this way we avoid 
situations in which the divergence force could become the 
main source of acceleration and thus instabilities due to tem- 
poral high numerical divergence, e.g. during strong interac- 
tions (see also section |4.3,2[) . This limitation of the diver- 
gence force also helps to maintain energy conservation (for 
details see Stasyszyn & Dolag (2011), in preparation). The 
basic method of divergence force subtraction is not conserva- 
tive, however, by limiting the correction we reduce possible 
transfer of energy associated with the numerical divergence 
to kinetic energy. 

Similarly to IKotarba et"afl l|2010r ), we apply radia- 
tive cooling, star formation and the associated supernova 
feedback, but exclude explicit supernova-driven galactic 
winds, f ollowing the sub-resolu t ion m ultiphase model devel- 
oped by ISpringel fc Hernquistl (I2003T), in which the ISM is 
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treated as a two-phase medi um ( McKee fc Ostrikerl 1 19771 
IJohansson fc Efstathioull2006t ). 

The implementation used in this paper has been tested 
in det a il (ISpringel et al.ll2qqil; ISpringellbOoBl ; ISpringel et al.l 
l2005bl ; IDolag fc Stasvszvnl l2009h . We refer the reader to 



total mass 


A/tot 


1.34 x 1O 12 M 


disk mass 


M diBk 


0.075 Mtot 


bulge mass 


A^bulgo 


0.025 JW tot 


mass of the gas disk 


A/gas 


0.2 Af disk 


exponential disk scale length 


Id 


8.44 kpc 


scale height of the disk 


I'D 


0.2 l D 


bulge scale length 


Id 


0.2 l D 


spin parameter 


X 


0.1 


virial velocity of the halo 


V v ir 


160 km s" 1 


half mass radius 


Rbnl{ 


«12 kpc 


half mass circular velocity 


"half 


^249 km s" 1 


half mass rotation period 


Thalf 


^295 Myr 


initial temperature 


Tdisk 


fslO 000 K 


initial magnetic field 


Bo 


- 10~ 6 G 


Multi-Phase model parameters 


gas consumption timescale 


tMP 


8.4 Gyr 


mass fraction of massive stars 


Pmp 


0.1 


evaporation parameter 


A 


4000 


effective SN temperature 


T SN 


4 x 10 8 K 


cold cloud temperature 




1000 K 



Table 1. Parameters of initial disk setup 



these studies for further details on the applied method. De- 
tailed discussions on the numerical divergence of the mag- 
netic field (V • B) in SPH simula tions of isolate d and inter- 
acting galaxies can be found in IKotarba et all l|2009l) and 
IKotarba et all (|2010m . 



3 SETUP 



3.1 Galaxies 



The initial conditions for the s piral galaxies a r e produ ced us- 
ing the me thod descr i bed b y ISpringel et alJ (|2005bh which 
is based on iHernauistl (|l993n . The galaxies consist of a cold 
dark matter halo, a rotationally supported exponential stel- 
lar disk, an exponential gas disk and a stellar bulge com- 
ponent. The pro files of the halo and bulge are based on 
IHernauistl (|l990h . The halo, stellar disk and bulge particles 
are collisionless N-body particles. The gas is represented by 
SPH particles. 

The parameters used for the initial setup of the three 
identical galaxies can be found in Table[T] The magnetic field 
in each disc is set to B x — Bo, disk and B y — B z = with the 
z-axis being the axis of rotation. The particle numbers and 
fixed gravitational softening lengths e for each galaxy can 
be found in Table [2] The minimum SPH smoothing length 
for the gas particles is thereby 7i§ph = 1-Oe = 25 pc/h, with 
h — 0.71 being the Hubble constant. This choice of param- 
eters results in particle masses of m g a S = widisk = Tibuige ^ 
1.19 ■ 10 5 M Q //i and m ha io « 21.43 • W 5 M Q /h, respectively, 
with the Hubble constant h = 0.71. More details on the 
properties, evolution an d stability of t h e disks evo lved in 
isolation can be found in IKotarba etakl (|2009l . 120101 ). 
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Component initial particle 


fixed gravitational 


number 


softening 


length e 




lP c ] 




Dark Matter 4.0 X 10 5 


110/h 




Disk - stars 4.8 x 10 5 


25/h 




Bulge - stars 2.0 X 10 


25/h 




Gas 1.2 x 10 5 


ne / J, 
ZD/ II 




Total 1.2 x 10 6 






(a) The Hubble constant is assumed to be h 


= 0.71 in this 


paper. 


Table 2. Particle numbers and softening 


f 1 en i a t n q 




initial separation [kpc//i] a 


' scp 


160 




r b = 

' scp 


320 


pericenter distance [kpc/h] a 


r a 


12 




p 


21 


Gl 


G2 


G3 


rotation around x (0) 45 


90 


30 


rotation around y (tp) 45 


DO 





rotation around z ((/>) 








initial v x [km s _1 ] -145 


73 


18 


initial v y [km s _1 ] -157 


-95 


63 



(a) The Hubble constant is assumed to be h = 0.71 in this paper. 

Table 3. Collision setup parameters 



3.2 Orbits 

We want to study a general case of a galactic interaction. 
Thus, the initial orbital setup of the three colliding galax- 
ies has been chosen arbitrarily without the aim of matching 
a particular observed system. In creating the initial condi- 
tions, the galaxies (including their magnetic field) are first 
rotated with respect to the plane of sky (ly-plane) around 
the x-, y- and z-axes by the angles 9, ip and (j), respectively. 
In order to guarantee a final merger of the three galaxies, 
two of them (galaxy 1 (Gl) and galaxy 2 (G2)) are set on 
a nearly parabolic Keplerian two-body orbit with an initial 
separation of the centers of mass of r^ ep and a pericenter dis- 
tance of Vp. Then, the third galaxy (G3) and the center of 
mass of the combined system Gl and G2 are set on a nearly 
parabolic orbit with r^ cp and r P . The values of the rotation 
angles, the initial separations, the pericenter distances and 
the initial velocities can be found in Table [3] 



3.3 IGM 

As an extension c ompared to our previous study in 
iKotarba et~aH (|201(]| '). we now also include an ambient IGM 
surrounding the galaxies in order to be able to realisti- 
cally study the magnetic field evolution in this IGM. The 
IGM is set up by placing additional gas particles in a 
regular hexagonal close-packed lattice (hep) arrangement. 
Close-packed lattice arrangements have been shown to be 
mor e stable for the par ticles than simple grid arrangements 
(see |Price~fc Bate 20071 and references therein). The particle 
mass is again ttiigm = m ga , B ~ 1.19 • 10 5 MQ/h. The vol- 
ume filled by the IGM is 500 x 700 x 300 (kpc/ft) 3 and 
centered at the center-of-mass of the initial setup of the 
discs. The density is pigm = 10 -29 g cm -3 (equivalent to 



~ 6 - 10 6 H-atoms cm 3 ), resulting in a particle number of 
iViGM = 43 x 70 x 31 = 93310. 

The initial morphology of our galaxies corresponds to 
fully evolved disk galaxies residing in their dark matter 
haloes, hence, we assume that the IGM in each scenario is 
already virialized. However, for simplicity, we assume a com- 
mon temperature for the IGM, which we set to the virial 
temperature at the virial radius of the simulated galactic 
haloes: 



{v 2 }fj,m p GM[im v 



3k E 



3r v i T ki 



Tigi 



6 • 10 b K, 



(1) 



where /i ~ 0.588 is the molecular weight (for fully ion- 
ized gas of primordial composition), m p the proton mass, 
fes the Boltzmann constant, G the gravitational constant, 
M the mass of the halo and r v i r its virial radius. We assume 
the initial magnetic field of the IGM gas to be directed in 
the x direction (where the x — y-plane is the orbital plane) 
with B x = Bo, igm- Thus, altogether, the galaxies are per- 
meated by a homogeneous magnetic field lying in the planes 
of the disks, and the IGM hosts a homogeneous field lying 
in the orbital plane. 

The addition of an ambient IGM has a further, numer- 
ical advantage: As the field is not dropping to zero at the 
disk edges, spurious calculations of the numerical divergence 
^sphV ■ B/|B| (where hsm is the SPH smoothing length) 
are avoided, resulting in higher numerical stability. 



3.4 Magnetic fields 

A detailed study of the cosmological seeding and evolution 
of magnetic fields in the early universe is beyond the scope 
of this paper. Rather, we want to gain better understanding 
of how magnetic fields might have evolved during the epoch 
of galaxy formation by considering three different scenar- 
ios of the magnetization of the colliding galaxies in the lo- 
cal universe, and, additionally, a further scenario excluding 
magnetic fields for comparison. 

The initial magnetic field strengths for each scenario 
are summarized in Table [4] With the first scenario (G6- 
IGM9), we aim to simulate a present-day galactic merger. 
Hence, within this scenario, we assume that the galaxies al- 
ready host an initial magnetic field of Bo, disk = 10 -6 G, 
and that the IGM is interspersed with an initial magnetic 
field of Bo. IGM = 10" 9 G. With our second scenario (G9- 
IGM9), we want to study the general situation of a com- 
mon magnetic field strength in the IGM and the galaxies. 
Within this scenario, the galaxies as well as the IGM host 
an initial magnetic field of 10 -9 G. In the third scenario 
(G9-IGM12), we assume that the magnetic field strengths 
are weaker by roughly three orders of magnitude than to- 
day, i.e. B ,disk = 10~ 9 G, and B ,igm = 10~ 12 G. In the 
last scenario (G0-IGM0), we exclude any magnetic field. For 
simplicity, all initial magnetic fields are assumed to be ho- 
mogeneous with only one non-vanishing component at the 
beginning of the simulations (which lies in the plane of the 
disks for the galaxies, and in the orbital plane for the IGM). 
The choice of the initial configuration is not important, as 
the timescales of the simulations are much longer than e.g. 
the turbulent timescales of the gas. Moreover, it takes more 
than 0.5 Gyr until the first encounter between the galaxies 
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Scenario 


-Bo, disk [G] 


#0,IGM [G] 


G6-IGM9 


lO" 6 


10" 9 


G9-IGM9 


10- 9 




G9-IGM12 


10- 9 


10- 12 


GO-IGMO 







Table 4. Initial uniform magnetic field strengths for the different 
scenarios in the plane of the disks (Bo, disk) an d m the ^-direction 
of the orbital plane (So,IGm)> respectively. 

and the magnetic field within the galaxies has enough time 
to redistribute and form a realistic configuration pr ior to the 
first encounter (see also lKotarba et al]|2009l . |2010| ). 

Finally, we let the system evolve for 200 Myr in order 
to allow possible numerical discontinuities introduced by the 
setup to relax before considering its physical properties. 



4 SIMULATIONS 

4.1 Morphological evolution 

Figs. [1] [2] and [3] show the mean line-of-sight total magnetic 
field |B|, the mean line-of-sight temperature T and the mean 
line-of-sight rms velocity v lnls , respectively, at nine different 
time steps for the G6-IGM9 scenario (i.e. the present-day 
mergerjj. To generate the images, the particle data in the 
500 3 (kpc/hf volume was binned to a 512 image using the 
code P-SMAC2 (Donnert et al., in p reparation), which ap- 
plies the gather approximation (see iDolag et al.l 120051 ). In 
each case, the color maps are overlaid with contours of the 
stellar surface density E* in order to indicate the stellar 
morphology of the system. The stellar density was binned 
usin g the TSC procedure (Tria ngular Shaped Cloud, see 
e.g. Hocknev fc Eastwood! Il988l ) . Here and hereafter, v vms 
of each particle is defined as the rms velocity around the 
mean velocity inside the SPH kernel of the particle (i.e. in- 
side its smoothing length /isph). This value is different from 
the rms velocity around the central velocity (i.e. the velocity 
of the particle considered) inside the kernel. The lat t er wa s 
defined as the "turbulent" velocity in iKotarba et"ai] < |201Ch . 
However, v Ilns is a more conservative estimator for the lo- 
cal turbulent velocity on the kernel scale. By averaging over 
all particles inside the kernel to obtain the mean we reduce 
the influence of SPH sampling noise, which is present on 
subkernel scales, on the estimate. Hence, we use « rms as the 
estimate of the turbulent velocity in this paper. 

At t = 300 Myr, the initial magnetic field has been 
redistributed and forms a realistic configuration in each 
galaxy (upper left panel in Fig. [T]). Due to the winding of 
the initially uniform magnetic field by differential rotation, 
the magnetic field in each galaxy forms a non-axisymmetric 
pattern, which can be recognized from the two magnetic 
arms extending from the galactic disc (best visible at the 
later times t = 654 M y r and t — 715 Myr in G3, see 
also iKotarbaet all 120091 . l20ld for a more detailed discus- 
sion of the winding process). As a consequence of their mu- 
tual gravitational attraction, Gl and G2 are moving towards 

1 See http://www.usm.uni-muenchcn.de/pcoplc/kotarba/public.html 
for the corresponding movies. 



each other and collide at t ~ 0.7 Gyr. Due to the lower 
ram pressure within the IGM than within the galaxies, the 
interaction-driven shocks (which are analyzed in more detail 
in sect ion [472]) are propagating favorably into the IGM, thus 
heating the IGM gas (upper central and right panels in Fig. 
0). Simultaneously, the rms velocity is enhanced and the 
magnetic field is strengthened in the shocked regions (upper 
right panels in Figs.[T]and[3]). During the collision, prominent 
tidal arms are developing (middle left panels in Figs. [T]- [3]). 
At t m 1.1 Gyr a second collision between Gl and G3 takes 
place, again accompanied by shocks and interaction-driven 
outflows. During the subsequent violent collision between 
G2 and G3 at t ~ 1.2 Gyr, further shocks and outflows are 
driven into the pre-shocked gas of the IGM (middle right 
panel in Fig. [3J|. The pre-shocked gas is now ejected by the 
cumulative shock fronts of this interaction, thus magnetizing 
the IGM. At t « 1.8 Gyr, the IGM within several 100 kpc 
around the collision debris is magnetized and highly turbu- 
lent (lower left panels in Figs. [T] and [3]). The debris continues 
to interact until their mutual final merger at t ~ 2.5 Gyr 
(lower right panels in Figs. Q]-[3]). By the end of the sim- 
ulation, the galactic magnetic field has an average value of 
approximately 10 -6 G, thus retaining its initial value, and 
the average IGM magnetic field reached a final strength of 
roughly 10 -8 G (lower right panel in Fig. [TJ. 

The global morphological evolution of the system within 
the G9-IGM9, G9-IGM12 and the GO-IGMO scenarios is sim- 
ilar to the evolution within the G6-IGM9 scenario presented 
above, i.e. the evolution of the stellar distribution does not 
change significantly, and the collisions and merger events 
take place at the same times0 For clarity, we do not show 
the corresponding figures to Figs. [1] - [3] for the other sce- 
narios here (see appendix EJ, but rather highlight the main 
differences below. 



4.2 Differences between the scenarios 

The most prominent differences between the different sce- 
narios are the evolution of the magnetic field and, simulta- 
neously, the behaviour of the shocks driven into the IGM 
during the interactions. Thereby, the shocks are propagat- 
ing faster (and gain higher Mach numbers, section l4.2.2|) for 
stronger initial magnetic fields. This is best visible during 
and shortly after the first collision (t « 700—900 Myr), when 
shocks are driven into the IGM for the first time. After the 
second collision, multiple shocks are propagating through 
the IGM which are additionally moving with different ve- 
locities depending on the scenario considered. Thus, even if 
it is possible to identify a single shock within one scenario, 
it is not possible to find the corresponding shock structure 
within the other scenarios. Therefore, within this section, we 
confine our analysis to the time before the second collision, 

2 This is not surprising, as it is commonly believed that the global 
morphological evolution of interacting galactic systems is deter- 
mined mainly by gravity and associated tidal action. However, it 
is interesting to note that in the early sixties, some authors be- 
lieved that the filaments and tails observed in interacting galactic 
systems have nothing to do with tidal phenomena and even con- 
sidered magnetic fields as the driving force for the morphologica l 
appearance of those systems j Vorontsov- Velvaminoll 196Cj , \l96%) . 
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Figure 1. The evolution of the magnetic field as a function of time for the G6-IGM9 scenario. Colours visualize the mean line-of-sight 
total magnetic field |B| (in units of 10 — 6 G) and contours correspond to stellar surface density E s tars- The contour levels are 10, 20 and 
50 M Q pc~ 2 . 
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Figure 2. The evolution of the temperature as a function of time for the G6-IGM9 scenario. Colours visualize the mean line-of-sight 
temperature T (in K) . Contours correspond to stellar surface density as in Fig. [T] 
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Figure 3. The evolution of the rms velocity as a function of time for the G6-IGM9 scenario. Colours visualize the mean line-of-sight 
rms velocity » Ims (in km s _1 ). Contours correspond to stellar surface density as in Fig. [T] 
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Figure 4. Shock propagation within the different scenarios, 
traced by the rms velocity. From top to bottom: the mean line- 
of-sight rms velocity at t = 715 Myr (left panels) and t = 859 
Myr (right panels) for the G6-IGM9, G9-IGM9, G9-IGM12 and 
the G0-IGM0 scenario, respectively. The color table is the same 
as in Fig. [3] whereby we only show the upper two galaxies Gl and 
G2. The higher the initial magnetic field (decreasing from top to 
bottom), the faster the shock propagation. The black boxes in- 
dicate the regions within which we analyse the shocks in more 
detail (see below). 



where we can identify and compare the shocks within every 
scenario. 



4-2.1 Differences in the shock propagation 

Fig- IU shows the shock propagation within the different sce- 
narios, which is best visibly traced by the rms velocity. From 
top to bottom: the mean line-of-sight rms velocity at t = 715 
Myr (left panels) and t = 859 Myr (right panels) for the G6- 



IGM9, G9-IGM9, G9-IGM12 and the G0-IGM0 scenario, re- 
spectively. The color table is the same as in Fig. [3] whereby 
we only show the upper two galaxies Gl and G2. 

The energy initially released in form of shocks by the 
subsequent interactions should be comparable within each 
scenario, as the setup of the galaxies and thus the impact 
parameters are the same. Nevertheless, Fig. [4] clearly shows 
that the stronger the initial magnetic field (decreasing from 
top to bottom), the faster the shocks propagate (note that 
the different scenarios differ only in the value of the initial 
field). The shock is the slowest for the scenario excluding 
magnetic fields (G0-IGM0, bottom panels). This behaviour 
suggests that the magnetic pressure associated with the 
magnetic field is able to additionally push the gas driven 
out by the interaction and thus accelerate the shocks. 



4-2.2 Differences in the magnetic field evolution 

Fig. [5] shows zoom-ins on the shock region at t = 715 Myr 
(black boxes in the left panels in Fig.[4| for the G6-IGM9 sce- 
nario (left panel), the G9-IGM9 scenario (central panel) and 
the G9-IGM12 scenario (right panel), respectively. Colours 
give the total magnetic field strength applying the same 
colour coding as in Fig. [T] Black vectors show the gas ve- 
locity whereby a length of ten on the given spatial scale 
corresponds to 500 km s _1 . The white line indicates the line- 
of-sight ri os used for the shock analysis presented below. The 
higher the initial magnetic field, the faster the shock (see also 
Fig. |4j , and the magnetic field is significantly strengthened 
behind the shock front. Independent of the initial magneti- 
zation model, the magnetic field strength behind the shock 
reaches values of the order of 10 -6 G, thus showing that the 
growth is more efficient for lower initial magnetic fields. 

Fig. [6] shows zoom-ins on the shock at t — 859 Myr 
(black boxes in the right panels in Fig- [4j , similar to Fig. [5] 
but with vector length of ten corresponding to 200 km s _1 . 
As before, the magnetic field strengths behind the shock 
reach a common value. However, within the tenuous IGM 
gas, the magnetic field reaches strengths of only 10~ 8 G. 

Fig. [7] shows the hydrodynamic values along the line- 
of-sight ri os for the shock at f = 715 Myr. The values 
are calculated using the SPH interpolation formalism. We 
show the shock properties for the G6-IGM9 scenario (blue 
lines), the G9-IGM9 scenario (red lines), the G9-IGM12 
scenario (black lines), and the G0-IGM0 scenario (green 
lines), respectively. From left to right and top to bottom: 
The number density n; the sound velocity c a , the Alfven 
velocity VAifVen — B/^/4-Kp and the magnetosonic veloc- 
ity c ms = a/c| + v Aif vell i respectively (whereby we do not 
show VAifven and Cms for the G0-IGM0 scenario); the gas 
velocity projected onto the line-of-sight v plo ; and, finally, 
the magnetic pressure P mag (except for the G0-IGM0 sce- 
nario), the thermodynamical pressure Ptherm and the shock 
energy density e s hock, respectively. We estimate the shock 
energy density by considering the difference between the 
projected downstream (ds) peak velocity Vp^ o peak (i.e. the 
shock velocity) and the average upstream (us) velocity Vp^ n , 

?i° r s n ) 2 . Magnetosonic Mach 
for the magnetized 



i.e. e shock = l/2p 
numbers M ms = 



C„ds 
V pro, peak 



(ll ds 
V pro, peak 



pro/ 

v us )/c u 



scenarios, and the Mach number M = (Upr ,peak 



for the G0-IGM0 scenario are given in the upper left panel. 
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Figure 5. Zoom-ins on the shock region at i = 715 Myr (black boxes in the left panels in Fig. [4| for the G6-IGM9 scenario (left), 
the G9-IGM9 scenario (central) and the G9-IGM12 scenario (right). Colours give the total magnetic field strength similar to Fig. [T] and 
black vectors show the gas velocity (a length of ten on the given spatial scale correspond to 500 km s" 1 ). The white line indicates the 
line-of-sight ri os used for the shock analysis. The magnetic field strength reaches f» 10 — 6 G behind the shock within each scenario. 




Figure 6. Zoom-ins on the shock region at t = 859 Myr (black boxes in the right panels in Fig. [3J|, similar to Fig. [5] but with a vector 
length of ten corresponding to 200 km s -1 . The IGM magnetic field strength behind the shock front reaches 10" 8 G within each 
scenario. 



Again, Fig. [7] clearly shows that the shock is propagat- 
ing faster the stronger the initial magnetic field (e.g. lower 
left panel). Thus, it is the fastest within the G6-IGM9 sce- 
nario, and the slowest within the GO-IGMO scenario. This 
increase in the shock intensity is also reflected in the pro- 
jected gas velocities (lower left panel) and the corresponding 
Mach numbers, which range from 1.7 for the GO-IGMO sce- 
nario to 6.0 for the G6-IGM9 scenariojf] Within each scenario 

3 Mach numbers of approximately 2-5 have also been found by 



(except for the GO-IGMO scenario), the magnetic field gets 
strengthened behind the shock (solid lines in the lower right 
panel). Thereby, the magnetic pressure behind the shock is 
of order of the energy density of the shock itself (between 
10~ 13 and 10~ 12 erg cm -3 ) within each scenario, thus con- 



Ijohansson et al. | ll2009ch in their numerical studies of gravita- 
tional heating through the release of potential energy from in- 
falling stellar clumps on galaxies. 
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Figure 7. Hydrodynamic values along the linc-of-sight ri os at t = 715 Myr for the G6-IGM9 scenario (blue lines), the G9-IGM9 scenario 
(red lines), the G9-IGM12 scenario (black lines), and the G0-IGM0 scenario (green lines), respectively. From left to right and top to 
bottom: The number density n; the sound velocity c s , the Alfven velocity ^Alfven an d the magnetosonic velocity c ms , respectively; the 
gas velocity projected onto the line-of-sight v pTO ; and, eventually, the magnetic pressure P mag , the thermodynamical pressure Pthcrm an d 
the shock energy density e s h oc k! respectively. Magnetosonic Mach numbers M ms for the magnetized scenarios are given in the upper left 
panel. For the GO-IGMO scenario, we give the Mach number M and do not show i>Aifvem c ms and P mag , respectively. 



firming that the growth of the magnetic field is the more 
efficient the lower the initial magnetic field. 

Fig. [8] shows the same quantities as in Fig. [Jj but for 
the shock at t = 859 Myr. Again, the shock is propagating 
faster the stronger the initial magnetic field. However, the 
Mach numbers are lower than at t = 715 Myr, ranging from 
1.4 for the GO-IGMO scenario to only 2.8 for the G6-IGM9 
scenario. Thus, also the shock energies are lower by roughly 
one order of magnitude (lower left panel). As before, the 
magnetic field gets strengthened behind the shock (lower left 
panel). However, for this weaker subsequent shock which is 
propagating through a pre-shocked gas further outside the 
galaxies, the growth is not as efficient as at t — 715 Myr. 
Hence, the magnetic pressure behind the shock does not 
reach the shock energy density. 

However, given the different sound speeds (or magne- 
tosonic velocities) within the different scenarios, also the hy- 
drodynamical timescales (~ c B /l, where I is a typical length- 
scale) are different. Thus, the comparison of the shocks at 
the same global dynamical time (defined by the gravitational 
interaction) might be questioned. In order to verify our con- 
clusions drawn on the basis of Figs. [7] and [8] we show how 
the shocks within each scenario compare when the density 
peaks of the shocks driven by the first collision have cov- 
ered the same distance (Fig. [9]). The corresponding different 
global dynamical times for each scenario are given in the 



lower left panel of Fig. [9] The lower the sound speed (upper 
right panel), the later in time the morphological agreement 
with the G6-IGM9 scenario (blue lines) is reached (upper left 
panel). The difference in the shock behaviour, i.e. the fact 
that the projected velocities are smaller for a smaller initial 
magnetization (lower left panel), is even more pronounced 
as before in Fig[7] This is reflected also in the Mach numbers 
(upper left panels in Fig. [7] and [5} and in the energies of the 
shocks (lower right panels in Fig.[7Jand|9]), and confirms our 
conclusion that the shocks are stronger for a stronger initial 
magnetization. 

In summary, when shocks and outflows are driven by the 
interactions, the magnetic fields get strengthened. Thereby, 
the weaker the initial magnetic field, the more efficient the 
growth of the field (except for scenario GO-IGMO, where we 
exclude magnetic fields). The magnetic fields are thereby 
most efficiently strengthened behind the shocks. By the end 
of the simulations, the magnetic field strengths and distri- 
butions within every scenario are comparable (see section 
14. 3. H and the appendix [AjJ. 

4-2.3 Differences in the temperature evolution 

According to the Rankine-Hugoniot shock jump conditions, 
the temperature behind the shock is proportional to the 
Mach number of the shock (more precisely, for 7 = 5/3, 
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Figure 8. Same as Fig. [7] but for the shock at t = 859 Myr. Note the different ^-ranges compared to Fig. [7] 



T ds /T us oc 5Af 2 + 14 - 3/M 2 , see lLandau fc Lifshit j 
1959). Thus, the higher Mach numbers within the different 
scenarios are reflected in the temperatures of the IGM 
behind the shock fronts (best visible for t — 715 Myr, 
Fig. I10p . Simultaneously, the region of shock heated gas 
is larger the stronger the shock (see also appendix IA2I) . 
Thereby, the IGM temperatures reach values up to 10 s 
K. Gas of this temperature can be expected to give rise 
to significant X-ray emi ssion by t h ermal bremsstrahlung. 
This was shown by e.g. ICox et all (|2006l ), who performed 
hydrodynamical simulations of mergers of gas-rich disk 
galaxies. They showed that the hot diffuse gas that is pro- 
duced by strong shocks attending the merger process can 
produce appreciable X-ray emission. Hence, it is interesting 
to note that the expected X-ray emission during galactic 
interactions increases with increasing magnetic field. 

As the nature of the presented galactic interaction is 
highly non-linear, our conclusions should be taken as gen- 
eral trends. The behaviour of the shocks in our simulations 
strongly indicates that in the presence of a magnetic field 
the shocks driven by an interaction are supported by the 
magnetic pressure, thus resulting in higher Mach numbers. 
Simultaneously, the magnetic field evolution in our simula- 
tions is a strong indication for a shock driven magnetic field 
amplification. As most of the gas in the Universe is mag- 
netized and shocks can be driven by different processes, we 
want to emphasize that the "magnetic shock acceleration" 
seen in our simulations might have implications within many 
astrophysical settings, e.g. galaxy clusters. 
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Figure 11. B rms = \J (B 2 ) as a function of time for the G6- 
IGM9 scenario (blue lines), the G9-IGM9 scenario (red lines), 
and the G9-IGM12 scenario (black lines). We separately show 
the IGM values (dotted lines) and the values inside the galaxies 
(solid lines). We distinguish between the IGM and the galaxies 
applying a density threshold of 10 -29 g cm -3 . The final values 
for the galactic and the IGM magnetic field, respectively, are the 
same within all scenarios. 
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Figure 9. Same as Fig. [7] but for these moments in time, at which the density peaks of the shocks driven by the first collision within 
each scenario have covered the same distance (the different times for each scenario are given in the lower left panel). The even greater 
difference in the Mach numbers compared to Fig. [7] confirms the conclusion that the shocks are supported by magnetic pressure. 
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Figure 10. Comparison of the temperatures at t = 715 Myr within the different scenarios. From left to right: the G6-IGM9, G9-IGM9, 
G9-IGM12 and the G0-IGM0 scenario. The color table is the same as in Fig. [2] The higher the Mach numbers (decreasing from left to 
right), the higher the temperature. 
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4.3 Global evolution 

4-3.1 Magnetic fields 

Fig. [11] shows the rms magnetic field B rms = 
\J (B| + By + _B|) as a function of time for the G6- 
IGM9 scenario (blue lines), the G9-IGM9 scenario (red 
lines), and the G9-IGM12 scenario (black lines). We sep- 
arately plot the IGM values (dotted lines) and the values 
inside the galaxies (solid lines). We distinguish between the 
IGM and the galaxies applying a density threshold of 1CP 29 
g cm -3 . 

Within all magnetized scenarios, the galactic magnetic 
field saturates at several p,G by the end of the simulation, 
independent of the initial magnetic field strength. The first 
collision between Gl and G2 initiates within all scenarios 
a steady growth of the galactic rms magnetic field, which 
continues throughout the second and third collision. This 
growth starts at a later time and is less efficient within 
the G6-IGM9 scenario (where the initial galactic field is 
already of order of the final value) than within the other 
scenarios. This indicates that it is easier to enhance a mag- 
netic field which is not yet of the order of the saturation 
value. In the subsequent evolution, the galactic rms mag- 
netic field strength within the G6-IGM9 scenario decreases 
again. Thus, at time of the intermediate merger (t « 1.8 
Gyr), the galactic magnetic field has a rms value between 
1 and 10 pG within all scenarios. The intermediate and the 
final mergers (t « 2 Gyr) do not lead to further growth of 
the field. 

The evolution of the IGM rms magnetic field (dotted 
lines) results in a common final value of approximately 1CP 8 
G within all scenarios. Again, the first collision initiates a 
growth of the IGM magnetic field. Within the G9-IGM12 
scenario, where the initial IGM field is four orders of mag- 
nitude smaller than the final value, this growth continues 
throughout the second and third collision until the final 
value is reached after approximately 500 Myrs. Within the 
G6-IGM9 scenario, the initial IGM field of 10~ 9 G grows 
by a factor of 10 at time of the first collision, thus reaching 
the final value shortly after this collision. None of the fur- 
ther collision and merger events lead to a significant further 
growth of the field. Most interestingly, the same initial IGM 
field of 10~ 9 G does not grow as efficiently at time of the 
first collision within the G9-IGM9 scenario. Rather, within 
the G9-IGM9 scenario, the IGM field grows more slowly and 
reaches the final value only after the final merger. The only 
difference between the G9-IGM9 and G6-IGM9 scenarios is 
the different initial galactic magnetic field, which is three or- 
der of magnitude higher within the G6-IGM0 scenario. Thus, 
the faster growth of the IGM field within the G6-IGM9 sce- 
nario shows that interaction-driven outflows transport mag- 
netic field energy from the galaxies into the IGM, resulting 
in a more efficient growth of the IGM magnetic field. 

The fundamental source of energy for the general 
strengthening of the magnetic field is the gravitational en- 
ergy released during the interaction of the galaxies. This 
gravitational energy is converted into kinetic energy of the 
particles, particularly, turbulence which is expected to be 
driven on kpc scales. The kinetic energy in turn is partly 
converted into magnetic energy by compression, shearing 
and folding of the magnetic field lines. This behaviour re- 
sults from the basic thermodynamical principle according to 



which the free energy of a system will always be distributed 
among all available energy channels. Given the high vari- 
ability of the simulated system and the limited resolution, 
however, detailed studies of the small-scale processes respon- 
sible for this distribution (e.g. the turbulent dynamo) are 
beyond the scope of this paper. Yet, simple magnetic field 
line compression would result in a tight correlation between 
the magnetic field strength and the gas density. Thus, it is 
possible to demonstrate the presence of the other processes, 
i.e. the shearing and folding of the magnetic field lines in 
the course of turbulent amplification, by analyzing the dis- 
tribution of the magnetic field strength as a function of the 
gas density. 

Fig. 1 121 shows B 2 /Bo as a function of p/po at three dif- 
ferent timesteps (from left to right: t = 20 Myr, t = 859 Myr 
and t = 2454 Myr) for the G9-IGM12 scenario (upper pan- 
els), the G9-IGM9 scenario (central panels) and the G6- 
IGM9 scenario (lower panels). The light blue dots corre- 
spond to particles with an initial density p ^ 10~ 29 g cm~ 3 
(i.e. galactic particles), and the residual particles (IGM) are 
marked in red. The "trimodality" of the distribution of the 
G9-IGM12 and the G6-IGM9 scenarios at t = 20 Myr is a 
result of the initial conditions: At the edges of the galaxies, 
IGM particles carrying a weak magnetic field overlap with 
galactic particles having a stronger field. Magnetic diffusion 
redistributes the field among these edge particles, leading to 
an increase of the IGM field and a decrease of the galactic 
field in these regions. This effect is of course not seen for the 
G9-IGM9 scenario, where all particles carry the same initial 
magnetic field. 

The black solid lines crossing all panels of Fig. 1 121 follow 
the equation B 2 / B 2 = (p/po) 4 ^ 3 , which is valid for isotropic 
compression of a turbulent magnetic field. If compression 
would be the only process responsible for the magnetic field 
growth, the distribution could evolve only along this line. 
Yet, the intermediate and final distributions of the mag- 
netic field strengths (middle and right panels) lie above this 
relation, showing that turbulent amplification has addition- 
ally strengthened the field (this is least prominent for the 
G6-IGM9 scenario, where the initial magnetic field strength 
is already of the order of the final strength). On the other 
hand, the final distributions (right panels) follow the 4/3 
inclination, thus showing that compression has to be a part 
of the overall strengthening process. The factor by which 
the IGM magnetic field (squared) gets strengthened is ap- 
proximately (10 4 ) 2 for the G9-IGM12 scenario and (10) 2 for 
the G9-IGM9 and the G6-IGM9 scenarios, respectively. The 
black dashed lines (right panels) correspond to the black 
solid lines but shifted by this factor (not fitted to the distri- 
bution), thus additionally demonstrating the importance of 
non- compressive processes for the magnetic field evolution. 

The horizontal solid and dashed dark blue lines in Fig. 
1121 (right panels) indicate the difference between, respec- 
tively, the initial galactic and IGM magnetic field (Bo, disk 
and -Bo, igm ) within each scenario, and the common final 
galactic rms value of 10~ 6 G (for the G9-IGM9 scenario, 
these lines lie on top of each other). As this final galactic 
rms value is the maximal magnetic field strength which can 
be acquired by both, the galactic and the IGM particles, 
the dark blue lines actually show the saturation levels of the 
magnetic field in the different scenarios. This is also visible 
from the "noses" on the right-hand side of the final distribu- 
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Figure 13. The mean (/iV ■ B/|B|) as a function of time for the 
G6-IGM9 scenario (blue lines), the G9-IGM9 scenario (red lines), 
and the G9-IGM12 scenario (black lines). We separately show the 
IGM values (dotted lines) and the values inside the galaxies (solid 
lines). We distinguish between the IGM and the galaxies applying 
a density threshold of 10 — 29 g cm -3 . The numerical divergence 
is small enough to guarantee numerical stability. 



tions, which show that particles which have already reached 
the saturation level can not acquire higher magnetic field 
strengths even if the density increases. 

In summary, Fig. [12] shows that the magnetic field 
growth driven by the interactions is due to both, the com- 
pression and the folding and stretching of the magnetic field 
lines in turbulent flows. It also confirms the saturation level 
of the magnetic field strength as shown in Fig. 1111 However, 
we note that the presented simulations are not designed to 
study the process of turbulent amplification in detail, where- 
fore our results should be interpreted in terms of general 
thermodynamics and magnetohydrodynamics. 



4-3.2 Numerical divergence 

Fig. [13] shows the mean numerical divergence measure 
(/isph|V ■ B|/|B|) as a function of time for the G6-IGM9 
scenario (blue lines), the G9-IGM9 scenario (red lines), and 
the G9-IGM12 scenario (black lines). We separately show 
the IGM values (dotted lines) and the values inside the 
galaxies (solid lines). A numerical divergence is expected 
to arise from the SPH divergence operator even for a field 
which is divergence free in the first place, but tangled on 
sub-resolution scales. This numerical divergence arises also 
in simulations where the magnetic field is expressed in terms 
of Euler potentials, which avoid physical divergence by def- 
inition. Comparisons of simulations using the Euler poten- 
tials with simulations using the direct implementation (as in 
this paper) have shown that the numerical divergence does 
not influence the evolution of the magnetic fie ld significantly 
as lo ng as the divergence measure is ^ 1 (jKotarba et al.l 
l2009t ). For the simulations presented here, the mean numer- 
ical divergence of all gas particles stays always below this 
tolerance value, except for the IGM value within the G9- 
IGM12 scenario during the third collision. Note, however, 
that this comparatively high divergence does not lead to 



an enhanced magnetic field growth (Fig. Generally, the 
higher the initial magnetic field, the lower the numerical 
divergence. This tre nd is consistent with previous studies 
jKotarba et alj|20icl ). The reason for this behaviour is the 
Lorentz force acting on the particles, which is always oppos- 
ing the motions leading to a change (growth) of the magnetic 
field. The stronger the field, the stronger this force. Thus, in 
the presence of a stronger field, random particle motions on 
sub-resolution scales are more efficiently suppressed. These 
motions result in a sub-resolution tangling of the magnetic 
field, and by reducing the motions, the numerical divergence 
is simultaneously reduced. Altogether, the numerical diver- 
gence should not influence the general conclusions on the 
magnetic field evolution presented in this paper (see also 
iKotarba etlIll2009l . [20Tol ). 

4-3.3 Pressures 

Saturation phenomena as seen for the magnetic fields in our 
simulations (Fig. Ill|) usually suggest some kind of energy 
equipartition. Given enough time, a thermodynamical 
system will always distribute its free energy to all the 
degrees of freedom available. A particular saturation value 
is thereby reached when the energy which is increasing 
balances the energy of the source responsible for the 
increase. For example, in case of the galactic dynamo, the 
main source for the amplification of the magnetic field is 
the energy of the particles which have a velocity component 
perpendicular to the galactic disc. The corresponding rising 
and descending particle motions are expected to become 
helical under the influence of the Coriolis for ce, resulting in 
the so-called a-effect (see e.g. iKulsrudl (|l999l ) for a review). 
Assuming these particles to be cosmic ray (CR) particles 
rises the expectation o f equipartition betwe en the magnetic 
and the CR pressure dHanasz et al.l l2009a[) . This e quipar- 
tition was recently shown by lHanasz et al.l (|2009bl ). More 
generally, in the framework of MHD, any motion of the gas 
leading to a growth of the magnetic field will be suppressed 
by the magnetic field itself via the Lorentz force as soon as 
the magnetic energy gets comparable to the kinetic energy 
of the gas. The magnetic energy is then converted into 
kinetic energy of the gas, thus maintaining equipartition 
between the magnetic and gas kinetic energy. In particular, 
the magnetic field is expected to be in equipartition with 
the turbulent energy of the g only velocity gradients 

can lead to a growth of the magnetic field. This concept 
of a "turbulent dynamo" is well know n from theory (see 
e.g. iBrandenburg fc Subraman ian 2005 for a review). Thus, 
within this section, we analyse our simulations with respect 
to the different pressure components within the galaxies 
and within the IGM. 

Fig. [14] shows the evolution of the volume weighted 
mean of different pressure components as a function of time 
for the G6-IGM9 scenario (upper panel), the G9-IGM9 sce- 
nario (central panel) and the G9-IGM12 scenario (lower 
panel). We plot the following pressures: magnetic pressure 
= B 2 /8tt (red lines), thermal pressure ftherm = l/2pc 2 
(green lines, with c s being the sound speed and p the gas 
density), turbulent pressure P tur b = l/2pv 2 ms (black lines) 
and ram pressure P ram = l/2pv 2 (blue lines, with v the 
total velocity of the gas particle considered, i.e. including 
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Figure 12. B 2 /B 2 as a function of p/po at three different timesteps (from left to right: t = 20 Myr, t = 859 Myr and t = 2454 Myr) for 
the G9-IGM12 scenario (upper panels), the G9-IGM9 scenario (central panels) and the G6-IGM9 scenario (lower panels). The light blue 
dots correspond to particles with an initial density p ^ 10 -29 g cm -3 (i.e. galactic particles), and the residual particles (IGM) are marked 
in red. The black solid lines follow the relation B 2 /B 2 , = (p/po) 4 / 3 ! which is valid for isotropic compression. The final distributions of 
the magnetic field strengths (right panels) lie above this relation, showing that compression is not the only process responsible for the 
growth of the magnetic field. The black dashed lines (right panels) show the same relation but shifted by the factor by which the IGM 
magnetic field has been strengthened within each scenario (not fitted to the distribution). The horizontal solid and dashed dark blue 
lines (right panels) correspond to the difference between, respectively, the initial galactic and IGM magnetic field (Bp. disk an d ^0 igm) 
within each scenario, and the common final galactic rms value of 10 — 6 G (for the G9-IGM9 scenario, these lines lie on top of each other). 
These are the saturation levels of the magnetic field in the different scenarios. 



turbulent components). Again, we distinguish between the 
IGM (dotted lines) and the galaxies (solid lines) applying a 
density threshold of 10 -29 g cm~ 3 . 

Within the galaxies (solid lines), the thermal and ram 
pressures within each scenario evolve smoothly with time, 
whereby the thermal pressure stays always below the ram 
pressure by roughly one order of magnitude. By the end of 
the simulations, these pressure components have the same 
values as at the beginning. The turbulent pressure (black 
lines) within each scenario increases by slightly more than 
one order of magnitude during the simulations, whereby each 
collision and merger event tends to increase the turbulent 
pressure. 

The galactic magnetic pressures (red solid lines) within 
each scenario evolve in a similar way as the turbulent pres- 
sures, increasing at the times of collision and merger events 
and partly decreasing afterwards (due to the dilatation of 



the shocked region and the corresponding dilution of the 
magnetic field energy). After the second collision, the mag- 
netic pressure is in approximate equipartition with the tur- 
bulent pressure within all scenarios. Before, the evolution of 
the magnetic pressures within the different scenarios differs 
because of the different initial magnetization. Within the 
G6-IGM9 scenario (upper panel) , where the initial magnetic 
pressure is already of the order of the turbulent pressure, it 
increases by roughly one order of magnitude during the first 
collision, thus even exceeding the turbulent pressure. How- 
ever, it decreases again to its initial value shortly after the 
first collision. Within the other scenarios, where the initial 
magnetic pressure is lower than the turbulent pressure, the 
galactic magnetic pressure increases during the first collision 
only up to the equipartition level, and does not decrease to 
its initial value after this collision. 

Within the IGM (dotted lines), the thermal and ram 
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Figure 14. Evolution of the volume weighted mean of differ- 
ent pressure components as a function of time for the G6-IGM9 
scenario (upper panel), the G9-IGM9 scenario (central panel) 
and the G9-IGM12 scenario (lower panel). We plot the follow- 
ing pressures: magnetic pressure P mag (red lines), thermal pres- 
sure Ptherm (green lines), turbulent pressure Pt ur b (black lines) 
and ram pressure P ram (blue lines). We distinguish between the 
IGM (dotted lines) and the galaxies (solid lines) applying a den- 
sity threshold of 10 -29 g cm -3 . By the end of the simulations, 
the magnetic pressure is of the same order of magnitude as the 
turbulent pressure within all scenarios. 



pressures within each scenario are near equipartition un- 
til the third collision. This is reasonable, as we assume the 
IGM temperature to be the same as the virial temperature of 
the haloes at the beginning of the simulations. However, at 
roughly the time of the third collision, the thermal pressure 
begins to decrease more efficiently than the ram pressure, 
and, in the subsequent evolution, stays below the ram pres- 



sure by a factor 10 - 30. This bahaviour is similar for each 
scenario. The IGM turbulent pressure (black dotted lines) 
differs form the IGM ram pressure by roughly two orders of 
magnitude during the whole simulation within each scenario, 
which is comparable to the ratios between the turbulent and 
ram pressures within the galaxies. 

The evolution of the IGM magnetic pressure (red dotted 
lines) follows the evolution of the turbulent pressure within 
each scenario after the second collision. This parallel evo- 
lution indicates a correlation between these pressure com- 
ponents. However, contrary to what is seen in the galaxies 
(solid lines), the IGM magnetic pressure stays by roughly 
one order of magnitude below the IGM turbulent pressure 
until the end of the simulations. This difference is most prob- 
ably a result of the numerical method by which we estimate 
the turbulence, within which shearing and inhomogeneous 
motions can lead to an overestimate of the turbulence. This 
can particularly happen within the IGM, where shocks prop- 
agate in a nongeneric direction. Thus, the similarity of the 
evolution of the IGM magnetic and turbulent pressures may 
be interpreted as equipartition between these components. 

Before the second collision, the evolution of the IGM 
magnetic pressure within the different scenarios again differs 
because of the different initial magnetization. Within the 
G6-IGM9 scenario (upper panel), the "quasi" equipartition 
with the turbulent pressure is reached already after the first 
collision, whereas it is reached only after the second collision 
within the other scenarios. Particularly, the IGM magnetic 
pressure within the G9-IGM9 (central panel) scenario does 
not grow as efficiently as within the G6-IGM9 scenario, 
although the initial IGM magnetization is the same within 
both scenarios. This difference is due to the stronger initial 
galactic magnetic field within the G6-IGM9 scenario: The 
strong galactic field is transported by interaction-driven 
outflows into the IGM, thus enhancing the IGM magnetic 
field more efficiently than within the G9-IGM9 scenario 
(see also section [4. 3. 111 . 

Summing up, we conclude that the saturation value 
within the galaxies of several /xG seen in all simulations 
corresponds to the equipartition between turbulent and 
magnetic pres s ure, a result which was also achieved by 
iKotarba et all l|2010t ). Also, the saturation value within the 
IGM (ss 10 -8 G) supports the assumption of equipartition 
between magnetic and turbulent pressure. 



5 SYNTHETIC RADIO EMISSION, 
POLARIZATION AND RM MAPS 

Radio observations of local galaxies and galaxy groups pro- 
vide maps of the distribution and structure of the magnetic 
fields in these systems. However, the radio maps are al- 
ways only one snapshot in time and observers face the prob- 
lem of explaining the origin of the observed magnetic field 
strengths and structures. For example, observations of the 
compact group of galaxies Stephan's Quintet, which consists 
of four interacting spiral galaxies and was recently modele d 
by means of JV-body simulations by iRenaud et all (|2010t ). 
show a prominent ridge of radio emi ssion crossing th rough 
the system in between the galaxies (|Xu et al.l 120031 ) . It is 
commonly believed that this ridge of radio emission corre- 
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sponds to a shock front driven by a former interaction be- 
tween two of the galaxies of the group. Due to the shock 
compression, the magnetic field of the ambient gas might 
have been amplified. Moreover, electrons are expected to 
get Fermi-accelerated within the shock, thus giving rise to 
the enhanced radio emission. 

It is the aim of this section to show how our simulated 
system of three merging galaxies might look like when ob- 
served at a typical radio frequency. As we are able to provide 
synthetic radio maps at every snapshot, we can assess the 
evolution and thus the origin of the magnetic field pattern 
seen in the synthetic maps. As detailed radio observations 
are available only for local groups of galaxies, we have cho- 
sen the G6-IGM9 scenario for the calculation of the synthetic 
maps. 

iKotarba et"afl l|201Ch have already presented such syn- 
thetic radio maps for an isolated galaxy and the Antennae 
system. Compared to their calculations, in this paper, we use 
a more precise formula for the synchrotron emission, which, 
however, does not change the results qualitatively (see e .g. 
Ginzburg fe Svrovatskiilll965l . iRvbicki fe Lightmadfl986l or 
Longairlll994l for more details on the standard theory of syn- 
chrotron radiation). Assuming an energy distribution of the 
relativistic CR electrons of a power-law form 

n{E)dE = nE~ p dE (2) 

(where n(E) is the number density of electrons, k a constant 
normalization factor and p the index of the power spectrum) , 
the synchrotron emissivity J v at a given frequency v is given 
by 



Jv 



m e c 2 (p + 1) 
V 4 12 



r["?-J- 

12 



(3) 



where F denotes the Gamma function. The magnetic field 
component perpendicular to the line of sight B± is taken 
from the simulation. The fre quency v and t h e ind ex p are 
input parameters. Similar to IKotarba et ail l)2010h we as- 
sume p = 2.6 and v = 4.86 x 10 9 H z, corresponding also to 
the values given bv lXu et all (|2003l ). The value of p gives a 
spectral index of the radio emission of a = (p — l)/2 = 0.1 



which is typic al for spiral gala xies jGioia et al. 19821) . an d 
given also by IXu et al] (|2003h and IChvzv fe Beckl (feooll ). 
Also, this value is close to the power-law slope of 2.7 of 
the all-particle spectrum of CRs in the energy rage of 
10 x 10 9 eV < E sC 3 x 10 15 eV (e.g. iBlasil 120081 and ref- 
erences therein). The constants are the speed of light c, the 
electron mass m e and the electron charge e. For a given CR 
energy density ecR, ft can be derived using Eqs. [2] 

e C R = / En(E)dE = k E x ~ v dE. (4) 

Similar to IKotarba et~al1 (|2010h we assume an energy range 
of -Emm = 10 9 eV to E 

max — 10 15 eV in our calculations. 
The cutoff at 10 9 eV is justified, since for a magnetic field 
strength of 10 [iG (which is the typical maximum field 
strength in our simulations) synchrotron emission is radi- 
ated mainly b y particles with an energy of « 7 x 10 9 eV 
(|Longairlll994l . eq. 18.51). For lower magnetic fields, the en- 
ergy of the electrons responsible for the observed radiation 
is even higher. 



CRs with energies between 10 9 eV and 10 15 eV are 
produced mainly through particle acceleration within the 
shocks of SN remnants. They propagate by advection and 
diffusion along the magnetic field lines throughout the 
galaxy or the IGM. Given the non-trivial ity of the physics 
involved in these processes (e.g. lBlasi|[200l) . the distribution 
of CRs within a galactic system is not easily modeled. How- 
ever, CRs are bound to the magnetic field, which, in turn, is 
coupled to the (ionized) thermal gas. Thus, as a reasonable 
but simple assumption, the CR energy distribution may be 
expected to be proportional to the thermal energy distribu- 
tion. Particularly, we assume that ecR = 0.01 ■ e t hcrm- 

J v is calculated for every particle, and the total syn- 
chrotron intensity /tot is subsequently obtained by integrat- 
ing along the line-of-sight. The polarized emission 7 po i, the 
predicted observed degree of polarization n o b s and the po- 
larization angles tp are calculated in the same way as in 
IKotarba et"al] (|201Ch . 

The calculations are implemented within the code P- 
Smac2 (Donnert et al., in preparation), which subsequently 
bins the values of /tot, / P oi, n o b s and ijiona grid using the 
gather approximation as before for Figs. |T| to [3] In order to 
account for the spatial isotropy of the emission, we multiply 
the calculated total and polarized intensities by the factor 



/obs 



4tt • d 2 ' 



(5) 



with d the assumed distance to the observer and d p i x the 
pixel size (corresponding to a beam diameter) of ~ 0.98 
kpc//i. We assume d — 80 Mpc, which is a typical distanc e 
for local groups, e.g. the Stephan's Quintet l)Xu et al.ll2003l ). 
Thus, dpi x corresponds to an angular resolution of « 1".3. 
For comparison, the maximum resolution of the Very Large 
Array (VLA) at v = 4.86 x 10 9 Hz (A « 6 cm) is 0".4, 
and for the Effelsberg telescope 2". 4. However, observations 
performed with these telescopes are usually presented with 
a lo wer resolution of « 5" — 15", e.g. the S tephan's Quin- 
tet (IXu et alj|2003h. the A ntennae galaxies (|Chvzv fc Beck! 
l2004h or M51 (|Beckll2009h . The upcoming Square Kilome- 



ter Array (SKA) is expected have 50 times the sensitivity 
of the VLA with a ma ximum angular resolution of 0".02 at 
1.4 GHz (A « 21 cm) (|Gaenslerll20"09h . 

Fig. [15] shows synthetic radio maps for the same nine 
time steps as in Figs. IA2I - IA6I Colours visualize the total 
intensity (in /ij/pixel). White contours show the polarized 
intensity, whereby the contour levels are 0.001, 0.01, 0.1, 1, 
10 and 50 /ij/pixel. The direction of the magnetic field is 
indicated by the black lines which are inclined by tp + n/2, 
whereby the data of the 512 2 array has been rebinned to 
a 50 2 array, thus lowering the assumed angular resolution 
to ~ 13". The length of this lines is scaled according to 
the degree of polarization n o b B , the length-scale is given in 
the lower right corner of each plot. We show the magnetic 
field lines only where the polarized intensity is higher than 
0.01 times the maximum polarized intensity in the beginning 
of the simulation, corresponding to a threshold of ~ 0.001 
/xjy/pixel. 

At the beginning of the simulation, the initially unidi- 
rectional magnetic field lines within the galaxies are wound 
up by the differential rotation of the disks. Hence, the mag- 
netic field vectors derived from polarization are showing a 
nearly toroidal pattern at t « 300 Myr (upper left panel). 
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Figure 15. Synthetic radio maps for the G6-IGM9 scenario at the same nine time steps as in Figs. IA2I - IA6I Colours visualize the 
total intensity (in /ij/pixel). White contours show the polarized intensity, whereby the contour levels are 0.001, 0.01, 0.1, 1, 10 and 50 
/ij/pixel. Magnetic field lines derived from calculations of polarization are shown in black. 
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After the first collision (upper central panel), when shocks 
are driven into the IGM, the total and polarized intensities 
are considerably enhanced behind the shocks (upper right 
panel). Thereby, the magnetic field exhibits a regular pat- 
tern aligned with the direction of the outflow. In the subse- 
quent evolution, prominent tidal arms are developing (mid- 
dle left panel). These tidal arms are also visible in the total 
and polarized radio emission and are traced by the mag- 
netic field lines. Presumably, the magnetic field lines have 
been stretched by the shear flows of the tidal structure. At 
time of and after the second collision (middle central and 
right panels), the polarized synchrotron emission is also vis- 
ible outside the disks, showing that the interaction driven 
shocks and outflows have already magnetized parts of the 
IGM. During the subsequent evolution (lower left panel), 
the gas motions inside the galaxies become more random 
and thus the polarization of the radio emission generally de- 
creases. At time of the final merger (lower central panel) 
shock driven gas flowing out of the merging system again 
gives rise to a high degree of polarization. By the end of 
the simulation, when most of the gas has been driven to the 
central core of the merged system (lower right panel), the 
total emission is concentrated around this core, and the po- 
larization is tracing the weak outflows which are still driven 
into the IGM. 

Fig. [15] does not reveal whether the polarized emission 
originates from an unidirectional or a reversing magnetic 
field configuration. Therefore, below we also present syn- 
thetic rotation measure (RM) maps of the simulated sys- 
tem. The value of RM gives the strength by which the po- 
larization vector of polarized radiation passing a magnetized 
plasma is rotated, whereby the rotation angle is cj> = RM-X 2 . 
The RM valu e for each simulated particl e is calculated ac- 
cording to (cf. iRvbicki fc Lightmanlll986T ) 

RM (particle) = -^^A n * B h ( 6 ) 

with n e the number density of thermal electrons (equal to 
the number density of thermal protons and thus propor- 
tional to the gas density) and Bu the magnetic field com- 
ponent along the line-of-sight. RM is positive (negative) for 
a magnetic field directed toward (away from) the observer. 
The cumulative RM is obtained by integrating along the 
line-of-sight, whereby individual RM values may add or can- 
cel. Thus, RM distributions which show reversals on small 
scales indicate a reversing magnetic field. The cumulative 
RM is binned on a grid as before in Fig. 1151 

Fig. [TB] shows RM maps for the G6-IGM9 scenario at 
the same nine time steps as in Fig. [15] Colors visualize the 
RM (in 10 -2 rad m -2 ) on a asinh-scale for a better visibil- 
ity of the small RM values in the IGM (caused by the low 
electron density in these regions). Note that on this scale 
values of I, 5 and 10 correspond to RM « 0.01, 0.74 and 
110 rad m -2 , respectively. 

In the very beginning of the simulation (t — Myr, not 
shown), the lower galaxy (G3), which is seen exactly face-on, 
does not show any RM because of the lack of a Bu compo- 
nent. The two other galaxies, which are inclined with respect 
to the line-of-sight, show a positive RM. The subsequent ini- 
tial infall of IGM gas onto the galaxies and the winding of 
the initial magnetic field due to differential rotation lead to 
the development of the RM patterns seen at t = 306 Myr 



(upper left panel). E.g., the RM is positive at both sides 
of the galaxy seen edge-on because the winding of the ini- 
tially homogeneous magnetic field results in a field pattern 
with the field lines directed towards the observer at both 
sides of the galaxy. The lack of RM in the central parts of 
the galaxies results from the cancelation of the (symmetric) 
contributions from both sides of the galaxies along the line- 
of-sight. Generally, the RM does not show reversals on small 
scales, showing that the magnetic field is unidirectional on 
scales ^ 10 kpc. This regularity is a result of the homo- 
geneity of the initial magnetic field. After the first collision 
(t = 654 Myr), the RM values within the regions of polar- 
ized emission (upper right panel, cf. Fig.ll5[) are all negative, 
strongly indicating an unidirectional magnetic field configu- 
ration. This is also true for the interaction-driven outflows 
at later times (middle left to lower right panels), whereby 
the sign of the RM may change depending on whether the 
outflow is driven predominantly towards or away from the 
observer. This large-scale homogeneity of the RM distri- 
bution is most probably a relic of the homogeneous initial 
IGM magnetic field. Furthermore, the RM tends to change 
sign on the edges of tidal arms (e.g. middle central panel, 
t — 1063 Myr), showing that the magnetic field within the 
arms is predominantly of galactic origin in contrast to the 
IGM magnetic field. At later times, the RM distribution 
within the galaxies and the tidal arms shows frequent rever- 
sals, originating in the reversing magnetic field tangled by 
interaction-driven turbulence. 

In summary, when shocks and outflows are driven out of 
the galaxies or tidal structures are forming, a high amount 
of polarized emission can be expected due to the stretching 
of the magnetic field lines by the gas flows. Thereby, the 
magnetic field structure may be unidirectional or reversing 
depending on the precedent magnetic field. Generally, high 
synchrotron intensity corresponds to high density regions, 
except for periods of intensive shock ejection. During these 
periods, the magnetic field is enhanced by the shocks and 
transported from the galaxies into the IGM by interaction- 
driven outflows. Q 

Given the importance of shocks for the synthetic polar- 
ization in our simulations it is reasonable to ascribe observed 
prominent regular magnetic field structures like the mag- 
netic field ridge observed in the Stephan's Quintet to shock 
activity. This assumption holds even better, as electrons are 
accelerated within shocks (which we do not model in our 
calculations), thus giving rise to an even higher synchrotron 
emissivity. 



6 CONCLUSIONS AND OUTLOOK 

We have presented for the first time high resolution simula- 
tions of a merger of three disk galaxies within an ambient, 
magnetized IGM. We have studied three different models for 
the initial magnetic field strengths within the galaxies and 
the IGM, respectively, and compared these models to a sim- 
ulation excluding magnetic fields. The initial magnetic field 
strength range form 10~ 12 G to 10 -9 G within the IGM, and 

4 Sec http:/ /www. usm.uni-mucnchcn.de/pcoplc/kotarba/public. html 
for a movie of the polarization and RM maps of the G6-IGM9 
scenario. 
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Figure 16. RM maps for the G6-IGM9 scenario at the same nine time steps as in Fig. 1151 Colors visualize the RM value (in 10 2 rad m 2 ) 
on a asinh-scale for a better visibility of the small RM values in the IGM. 
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from 1CP J G to 10 G within the galaxies. We find that 
the magnetic field saturates at a value of several fiG within 
the galaxies and at roughly 10~ 8 G within the 1GM, inde- 
pendent of the initial magnetic field. This saturation levels 
correspond to equipartition between the magnetic and the 
turbulent pressure in t he system. This result is in agreement 
with previous studies of lKotarba et al.l (|2010i ) , who have pre- 
sented simulations of the interaction of two disk galaxies, 
partic ularly the Antennae system. However, iKotarba et al.l 
(|2010l ) did not include an ambient IGM, wherefore they 
could not study the behaviour of shock propagation within 
the IGM and its magnetization. The simulations presented 
in this paper show that the shock propagation within the 
IGM is changed significantly depending on the initial mag- 
netic field model. Thereby, the stronger the initial magnetic 
field, the stronger the shocks driven into the IGM. This re- 
sult suggests that the shocks are supported by the magnetic 
pressure, resulting in higher Mach numbers in the presence 
of a strong magnetic field. 

The main findings presented in this paper can be sum- 
marized as follows: 

• The magnetic field within the galaxies grows to the 
equipartition value between the turbulent and magnetic en- 
ergy density during the subsequent interactions between 
the merging galaxies. The IGM is magnetized by outflows 
and multiple shocks driven by the interactions up to nearly 
equipartition with the turbulent energy density within the 
IGM. 

• The final saturation value for the galaxies is of order 
of several /xG, and the final IGM magnetic field strength 
has an average value of 10~ 8 G, independent of the applied 
initial magnetic field. As the setup of the system was chosen 
arbitrarily and the system was allowed to evolve freely in 
time, it is very interesting that the final magnetic field values 
are consistent with observations. 

• The growth of the IGM field is more efficient, the higher 
the galactic magnetic field, suggesting that magnetic energy 
is transported from the galaxies into the IGM. 

• The initial values of the magnetic field in the galaxies 
and in the IGM affect the propagation of shocks in the IGM: 
The stronger the initial field, the faster the shock propa- 
gation, suggesting that the shocks are gaining higher Mach 
numbers due to magnetic pressure support. This effect might 
be referred to as "magnetic shock acceleration" . 

• The higher Mach numbers are also reflected in the 
higher temperatures of the shock heated IGM gas, and the 
shock-heated region is larger the stronger the shock. 

• Shocks play an important role for the polarized emission 
of an interacting system. Always when shocks are driven 
into the ambient IGM by an interaction, a high amount of 
polarized emission can be expected. 

The presented simulations of a merger of three disk 
galaxies do not only provide insights in the evolution and 
significance of magnetic fields in a highly nonlinear envi- 
ronment, but also agree well with observations. The typi- 
cal observed value of several fiG in different types of local 
galaxies shows to be a typical value which arises naturally in 
merging systems. Also, the typical observational estimates 
for the upper li mit of IGM magnetic fields of approximately 
10" 9 -10" 8 G (jKronberg et alj|2008l and references therein) 



are in agreement with the IGM magnetic field in our simu- 
lations. 

Given the complexity of the simulated system, for a 
future project, it would be worthwhile to perform detailed 
shock simulations in a magnetized IGM setting. Parameter 
studies of the dependence of shock characteristics on the 
applied magnetic field strength and structure would lead to 
a deeper understanding of the evolution of magnetic fields 
during galactic interactions. 

Finally, we emphasize that the efficient strengthening 
of magnetic fields during subsequent galaxy interactions up 
to levels consistent with observations is of particular inter- 
est within the framework of CDM clustering models. On 
th e basis of the simula tions presented in this paper and 
in IKotarba et alj (|2010T ) , we can assume that every galac- 
tic interaction contributes to the magnetization of the af- 
fected galaxies and the ambient IGM. As the phase of struc- 
ture formation in the Universe is accompanied by frequent, 
subsequent interactions of galaxies and galactic subunits, 
the observed present-day magnetic fields might at least to 
some extent be the result of interaction-driven amplifica- 
tion processes in the early universe. Furthermore, the trans- 
port of magnetic energy from the galaxies into the IGM 
by interaction-driven outflows could explain the existence of 
IGM magnetic fields of order 10~ 9 - 10~ 8 G already at high 
redshifts. We note, however, that the simulations presented 
here picture the interaction of three fully evolved present- 
day galaxies, wherefore they can not be compared directly 
with interactions in the early universe. High redshift galax- 
ies are known to be very different than present-day galaxies. 
However, subsequent galaxy interactions were more frequent 
in the early universe. Thus, it is reasonable to assume that 
the cumulative interaction-driven magnetic field amplifica- 
tion of even small initial seed fields might have magnetized 
young galaxies and their environment already at high red- 
shifts. More studies of magnetic field evolution during the 
early phases of the Universe, preferably in a cosmological 
context, would thus be of particular interest for the whole 
astronomical community. 
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APPENDIX A: APPENDIX 

In Figs. [1] [2] and [3] we showed the evolution of the magnetic 
field, the temperature and the rms velocity as a function of 
time for the G6-IGM9 scenario. Below, we show the corre- 
sponding figures for the G9-IGM9, G9-IGM12 and G0-IGM0 
scenarios, respectively (see table [3J. 

Al Magnetic fields 

Figs. IA1I and IA2I show the evolution of the mean line-of- 
sight magnetic field as a function of time similar to Fig. [TJ 
but for the G9-IGM9 and G9-IGM12 scenarios, respectively. 
The final distributions and strengths of the magnetic fields 
are comparable within every scenario. 

A2 Temperatures 

Figs. IA3I IA4I and IA5I show the evolution of the mean line- 
of-sight temperature as a function of time similar to Fig. [1] 
but for the G9-IGM9, G9-IGM12 and G0-IGM0 scenarios, 
respectively. The stronger the shocks driven into the IGM, 
the higher the temperatures. 

A3 RMS velocities 

Figs. IA6I IA7I and IA5I show the evolution of the mean line- 
of-sight rms velocity as a function of time similar to Fig. [TJ 
but for the G9-IGM9, G9-IGM12 and G0-IGM0 scenarios, 
respectively. The higher the initial magnetic field, the faster 
the shock propagation within the IGM. 



Magnetic fields in colliding galaxies 



log( M [10* G] ) 
-2 -1 







I I I I I I I I | I I I I I I I I I | I I I I I I I I I | I I I I I I I I I 



I I I I I I I I I | I I I I I I I I I | I I I I I I I I I | I I I I I I I I I 



I I I I I I I I | I I I I I I I I I | I I I I I I I I I | I I I I I I I I 



300 



100 



t = 30S Myr 



Gl 



100 



-300 



200 



100 



G2 



G3 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



t = 654 Myr 



Gl 



G2 



G3 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



t = 715 Myr 
Gl 



G2 




G3 



i n ii i| i n ii ii n 1 1 ii i ii n i 



859 Myr 



Gl 



t = 1063 Myr 



t = 1306 Myr 



G2 



100 



-300 



200 



100 




G3 




G3 



G3 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
t = 1779 Myr 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
t = 3045 Myr 



i ii ii n ii | ii n n ii i| i n ii ii n 1 1 ii i ii n i 
. t = 3454 Myr 






100 1 



-300 



Gl 



■ ' I I ' 



'I I ■ ' 



'I I I 



-100 100 
x [kpc] 



-100 100 
x [kpc] 



100 100 
x [kpc] 



Figure Al. Same as Fig.Q] but for the G9-IGM9 scenario. 
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Figure A2. Same as Fig.Q] but for the G9-IGM12 scenario. 
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Figure A3. Same as Fig. E] but for the G9-IGM9 scenario. 
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Figure A4. Same as Fig. [2] but for the G9-IGM12 scenario. 
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Figure A5. Same as Fig. [2] but for the G0-IGM0 scenario. 
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Figure A6. Same as Fig.E] but for the G9-IGM9 scenario. 



Magnetic fields in colliding galaxies 




Figure A7. Same as Fig. [3] but for the G9-IGM12 scenario. 
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Figure A8. Same as Fig.|3] but for the GO-IGMO scenario. 



